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We present measurements of the superconducting critical temperature T c and upper critical field 
H C 2 as a function of pressure in the transition metal dichalcogenide 2H-NbS2 up to 20 GPa. We 
observe that T c increases smoothly from 6 K at ambient pressure to about 8.9 K at 20 GPa. 
This range of increase is comparable to the one found previously in 2H-NbSe2. The temperature 
dependence of the upper critical field H C 2(T) of 2H-NbS2 varies considerably when increasing the 
pressure. At low pressures, H C 2(0) decreases, and at higher pressures both T c and -ff C 2(0) increase 
simultaneously. This points out that there are pressure induced changes of the Fermi surface, which 
we analyze in terms of a simplified two band approach. 

PACS numbers: 74.70.Xa, 74.62. Fj, 74.70.Xa 
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INTRODUCTION 

2H-NbSe2 belongs to the family of transition metal 
dichalcogenide compounds and presents a charge density 
wave (CDW) below Tqdw — 33 K, which coexists with 
superconductivity (T c = 7.2 K) 04I3 ■ 2H-NbS 2 is a 
related two-band superconductor, with similar T c , and 
no charge order [ll|, ll2j . The crystal structure of these 
materials consists of hexagonal transition metal - chalco- 
gen sandwiches which are coupled through weak van der 
Waals forces, leading to hexagonal layers with large c- 
axis constant and strongly anisotropic properties. Com- 
pressibility is larger along the c-axis than in- plane [131 ] - 
Changes in the electronic properties are produced by al- 
tering the lattice constants, using compositional tuning 
(substitution 14 1 irradiation [la , Il6| or intercalation be- 
tween layers 17Hl9j|) and applying pressure. Pressure 
has been shown to lead to an increase of the critical tem- 
perature in 2H-NbSe2 with a maximum T c around 8.5 K 
at 10 GPa. CDW critical temperature decreases to zero 
at 5 GPa in this material, and shows significant pressure 
induced modifications in other compounds of the same 
family [13j, I2CH24 1 . The importance of local strain has 
been highlighted recently [25j. In Ref.[20j, authors apply 
pressure to 2H-NbSe2 and find that the effective dimen- 
sionality of the electronic structure is increased above 4.6 
GPa. 

To characterize the electronic changes suffered under- 
pressure in bulk materials, the measurement of the up- 




FIG. 1: Structure of 2H-NbS2- The niobium atoms are surrounded 
by sulphur atoms following a trigonal prismatic coordination. The 
S-Nb-S layers that make up each packet are covalently bound. The 
coupling between pairs of packets is dominated by van der Waals 
forces. The lattice parameters are a = 3.321A, b = 5.751A and 
c= 11.761A [H. 



per critical field, H C 2 is a simple and useful tool. In clean 
superconductors (with a mean free path, £, greater than 
the coherence length, £), the usual dependence of H c2 
on temperature is given by Helfand-Werthammer theory, 
which assumes a simple single band spherical Fermi sur- 
face [27j. It consists of a linear increase of H C 2 close 
to T c , which flattens out at zero temperature. Within 



this theory, H c2 (0) 
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quantum, 7 w 0.577 Euler's constant and vf the Fermi 
velocity. The slope of the linear increase of H C 2 close to 
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detailed treatment for complex Fermi surfaces developed 
in Refs.[30l l31| . shows that the upper critical field can 
have a strong dependence on temperature, which allows 
determining the Fermi velocity vp and electron-phonon 
coupling A parameters on different parts of the Fermi 
surface. 

In the two-gap compounds MgB2, YM2B2C and 2H- 
NbSe2, it was found that the ambient pressure H C 2(T) 
has a strong positive curvature close to T c . The asso- 
ciated difference in the Fermi velocities in both bands 
leads to differing coherence lengths for each band and 
thus two characteristic features in the upper critical field, 
giv ing t he p ositive curvature observed in the experiment 
[23j, [24j, |30j, [32J. In these compounds, the form of such 
positively curved H c2 (T) has a strong dependence as a 
function of pressure, from which electron phonon cou- 
pling and Fermi surface velocities have been obtained. 
In 2H-NbS2 , scanning-tunneling microscopy also revealed 



the existence of two superconducting gaps [Tl|. Sub- 
sequent heat capacity measurements also show two-gap 
superconductivity and a small positive curvature of the 
upper critical field [12j. Moreover, the temperature de- 
pendence of the superfluid density measured by lower 
critical field and magnetic penetration depth is very sim- 
ilar to 2H-NbSe2, and is well described by a two gap 
model[33l.l34j|. Here, we present new measurements of T c 
and H C 2 (T) as a function of the applied pressure in 2H- 
NbS2- We observe a smooth increase in T c as a function 
of pressure. H C 2(T) has a positive curvature at ambient 
pressure, with an anomalous pressure dependence, which 
evidences pressure induced changes in the Fermi surface. 



Its diameter and depth allowed the insertion of the sam- 
ple as well as of some ruby balls that were used to deter- 
mine the pressure through the ruby fluorescence method 
[361 ] . Pressure was transmitted by a methanol-ethanol 
mixture (4 : 1), which has given quasi-hydrostatic con- 
ditions up to the pressures of interest in our experiment 
23l l24j . |37| . The susceptometer is the same already de- 
scribed elsewhere [23(. It was designed to obtain the 
largest signal to-noise ratio, with a pickup coil wounded 
very close to the sample space. Two primary-secondary 
coils systems of about 4 mm in diameter and 2 mm in 
height were used. The first one was located surrounding 
the sample, which was then within the field created by 
the primary, whereas the secondary acquired the voltage 
due to any change occurring in sample's properties. The 
other primary-secondary assembly was glued beside, with 
no sample inside. The two primaries were connected in 
series, so the resulting magnetic field was the same in 
both cases. The secondary coils were connected in series 
opposition, so the large signals due to the secondary coils 
themselves could be removed from the beginning. After 
further compensation by means of an attenuator and a 
phase shifter the total signal was detected with a lock-in 
amplifier. For each applied pressure, T c and H C 2 were 
obtained by measuring the magnetic susceptibility as a 
function of temperature and at different magnetic fields, 
applied parallel to the c-axis. The critical temperature 
and field were determined by the onset of the supercon- 
ducting transition curves, defined as the intersection of 
two tangents, one to the flat portion of the curve above 
and the second to the steepest variation in the signal 
below the superconducting transition. 



RESULTS 



EXPERIMENT 

We have measured single crystalline samples of 2H- 
NbS2 , with T c around 6 K and a residual resistivity ratio 
around 10. They were grown as described in Ref.[35j, 
and have a hexagonal shape, with lateral dimensions of 
about 150 /im and thickness around 30 fim. We loaded 
them into a pressure cell made from copper-beryllium al- 
loy, within the sample space delimited by the diamond 
anvils, which are 0.7 mm in culet diameter, and a gas- 
ket made from NiMo alloy. The diamond anvils were 
mounted onto a couple of sapphire cylinders inserted into 
the bores of two Cu-Be pieces. This choice of materials 
guarantees that the inductive coupling between the coils 
and the neighboring parts of the cell can be taken as neg- 
ligible. The gasket had an initial thickness of 300 /j,m, 
which was later reduced to about 60 /im in its center, af- 
ter indentation. A small orifice was made by means of arc 
discharges between the gasket and a molybdenum needle. 



We find an ambient pressure superconducting critical 
temperature of 5.7 K. In previous measurements, no no- 
ticeable increase was measured below 1 GPa 38I440J . Fig. 
E^a) shows the variation with temperature of the suscep- 
tibility in 2H-NbS2 sample, for applied pressures ranging 
between GPa and 19.8 GPa. Fig. Mb) displays the 
variation of T c as a function of pressure. Clearly, there is 
a progressive increment of T c with pressure. A T c maxi- 
mum is likely to exist, but above 20 GPa. Below 9 GPa 
T c increases with a slope dT c /dP — 0.09 K/GPa, which 
further grows to 0.22 K/GPa between 9 and 14, and then 
decreases to 0.16 K/GPa for higher pressures. 

The magnetic field dependence of the susceptibility un- 
der pressure is shown in Fig[3]for different temperatures 
and at 3 GPa. There is a smooth evolution of the sus- 
ceptibility with magnetic field, from which we can easily 
extract the upper critical field as the onset of the tran- 
sition. The transition widens significantly at lower tem- 
peratures, as expected for type II superconductor. The 
form of H C 2(T), discussed in the following figure, does 
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FIG. 2: (a) Variation of the susceptibility of a 2H-NbS2 sample as 
a function of temperature and for several applied pressures. The 
arrow and the lines show the way we used to extract the corre- 
sponding value of T c . (b) Critical temperature as a function of the 
applied pressure. Triangles and full circles correspond to increas- 
ing and decreasing pressure, respectively. Lines are a guide to the 
eye and simply join data points. Saturation is eventually expected 
above 20 GPa. 



not depend on chosing H C 2(T) from onset, midpoint or 
lower part of the transition. 

Previous ambient pressure results on the upper criti- 
cal field of 2H-NbS2 yielded a positive curvature above 
0.5T C , and a zero-temperature critical field value of about 
2.6 — 2.7 T 12|, |4l|. Our measurements of the ambient 



pressure temperature dependence of the upper critical 
field of 2H-NbS2 are shown in the upper panel of FigHJ 
and essentially confirm previous findings. We observe a 
slight positive curvature with an extrapolated zero tem- 
perature critical field of about 2.6 T. This positive curva- 
ture, although less pronounced, is similar to that found 
in H c2 {T) in 2H-NbSe2- When applying pressure in 2H- 
NbS2, the critical temperature increases, but the upper 
critical field at low temperatures decreases strongly up 
to 8.7 GPa, above which it increases together with T c . 
Our results show that H C 2(0) drops by a factor 1.5 be- 
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FIG. 3: Variation of the susceptibility of a 2H-NbS2 sample as 
a function of the magnetic field, for a pressure of 3 GPa and for 
several temperatures. The arrow and the lines show the way we 
used to extract the corresponding value of H C 2(T). Note that the 
transition width significantly increases with magnetic field. Other 
criteria for determining the transition temperature, such as mid- 
point or lower part of the transition region, lead to similar H C 2(T) 



tween and 8.7 GPa and then rises by a factor of 1.5 
for highest attained pressures around 20 GPa. This is 
a peculiar behavior. If Fermi surface parameters do not 
change, theory predicts that H c2 (0) oc T c 2 [13, [H] . 

Thus, the non-monotonous variation of H C 2(T) is at 
odds with most simple single band BCS theory, and we 
explored how a two band scenario could explain this 
behaviour. For that purpose, we used the same type 
of calculations of the upper critical field of strong cou- 
pling (multiband) su per conductors and its pressure de- 



pendence as in [23, [24|, |42|, |43[, and numerically lin- 
earize H C 2(T) equations obtained from Eliashberg the- 
ory 4J] . The main change of H C 2 (T) under pressure oc- 



curs in the form of the positive curvature and the value 
of H C 2 (0) , which can be used to determine values of the 
Fermi surface properties, by parametrizing the Fermi sur- 
face in two main different subgroups of electronic exci- 
tations. We consider two subgroups of electrons, with 
bare Fermi velocities vfi and i>F2, and coupling param- 
eters Ay (i,j = 1,2, with subindex 1 for the electronic 
group with largest pairing strength) . There is of course 
no possibility to deduce a unique set of parameters sim- 
ply from our H C 2(T, P) data: we rather choose to find the 
simplest possible scenario which fits the whole pressure 
dependence of H C 2(T). 

The choice for our scenario was guided by the above 
mentioned peculiar behavior, namely that T c increases 
with pressure, whereas the slope of H C 2(T) close to T c 
decreases. An increasing T c suggests an increasing par- 
ing strength, but at the same time, the decreasing slope 
suggests an increased Fermi velocity, and thus, a decreas- 
ing effective mass, which is unexpected with increased 
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FIG. 4: Temperature dependence of the upper critical field H C 2(T) 
as a function of pressure (empty and filled symbols). Fits to the 
model explained in the text are shown as lines. Note the decrease 
in the upper critical field until 8.7 GPa (upper panel), which ceases 
above this pressure (lower panel). 



pairing interactions, or an increased Fermi wave vector. 
In the latter case, we would also expect an increased 
bare density of states, for instance through an increased 
Fermi pocket volume. So we checked if a simple sce- 
nario, where the mere increase of the Fermi surface vol- 
ume of the main electronic group, at the expense of a 
decreased Fermi surface volume of the second electronic 
group, could be suitable, with a larger Fermi surface 
(smaller slope of -ff c2 ) for the first group to reproduce 
the change of slope with pressure. To be quantitative, 
we introduce two (and only two) parameters to describe 
the pressure evolution of both T c and H c2 (T) : these pa- 
rameters are pi{P) and P2{P), which can be thought 
as the ratio of the bare density of states of each elec- 
tronic group under pressure with respect to the density 
of states at zero pressure, due to the respective Fermi 
pocket volume change. We can then deduce the pressure 
evolution of the electron-phonon coupling parameters as: 
Aij(-P) = ^ij(®)Pj(P)> an d f° r the bare Fermi velocities: 



vfj(P) = vpi(0)pi(P). The renormalized Fermi veloc- 
ities Ejl are as usual: v* Fi (P) = v Fl (P) — ^ 4 . We 

fixed the values of the mean phonon frequency, 9 = 54.5 
K and Coulomb pseudopotential, p* = 0.1, and assume 
that they are pressure independent. We find that indeed, 
as shown in FigUl the complete evolution of H C 2(T,P) 
can be fitted when starting from the initial set of values: 
Au(0) = 1, A 22 (0) = 0, A 12 (0) = 1.1, A 21 (0) = 0.55, 
v F1 (0) = 3.1 10 5 m/s (v Fl (0) = 10 5 m/s), v F2 (0) = 0.155 
10 5 m/s (f^ 2 (0) = 10 4 m/s). With these initial values of 
the electron phonon coupling constants, pairing is con- 
trolled by the first electron group and by the interac- 
tion between the first and second electron groups. Fig[5K 
shows the pressure evolution of the renormalized Fermi 
velocities, and FiglSb, that of the coupling parameters 
or equivalently, of p\ and p 2 . This demonstrates that 
H C 2(T, P) (including T C (P)) is compatible with a scenario 
where one of the Fermi pockets expands while the other 
shrinks, and with a constant pairing potential. Only the 
bare density of states related to the Fermi pocket volume 
is changing with pressure. 



DISCUSSION AND CONCLUSIONS 

Our results show that the variation of the Fermi surface 
parameters is monotonic (up to 20GPa), and there are no 
maxima nor minima as in 2H-NbSe 2 
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dome shape of T C (P) is found which peaks around 10 GPa 
and the CDW disappears at 5 GPa [20j , where a kink in 
T C (P) is found. The latter is associated to changes in one 
of the Fermi surfaces, with a possible Lifshitz transition 
due to one band shifting below the Fermi level. In 2H- 
NbS 2 , the strong decrease of the zero temperature upper 
critical field below 8.7 GPa, associated with a slight in- 
crease in critical temperature, can be explained by a de- 
crease in the renormalized Fermi velocity and an increase 
in the electron-phonon coupling of one part of the Fermi 
surface (see Fig. [5]), due to Fermi surface modifications. 
This produces a more pronounced positive curvature in 
H C 2(T), as the differencies in Fermi surface parameters 
increase. Above about 10 GPa, variations are smoother 
and essentially governed by slightly increasing renormal- 
ized Fermi velocities of the rest of the Fermi surface. 

The Fermi surface of 2H-NbS 2 is not known in detail, 
but has likely the features which are believed to be com- 
mon to similar transition metal dichalcogenides, namely 
two pairs of concentric cylindrical sheets derived from Nb 
4d electrons 45j,|46|. Fig. [5] shows that pressure induced 
modifications seem to saturate near 20 GPa, therefore 
indicating that T c cannot be expected to be much higher 
than 9 K. It also hints to a shrinking of a part of the 
Fermi surface, while the other one grows at its expense, 
as vf i evidences. These features of the superconduct- 
ing properties cannot be related to a competition with a 




FIG. 5: Pressure dependence of the renormalized Fermi surface 
velocities v%,.(P)/vp.(0) (a) and of the electron-phonon coupling 
parameters Xij (P) / \ij (0) = pj{P) (b) for both subgroups of elec- 
trons in 2H-NbS2- Lines simply join points and are a guide to the 
eyes. 



CDW, as charge order is not present in 2H-NbS2. It has 
been suggested that the absence of CDW order in this 
compound is due to anharmonic effects, and that the su- 
perconducting properties are essentially determined by 
the anisotropy and strength of the electron-phonon cou- 
pling |47J |. 

2H-NbS2 is at the verge of CDW, which is favored by 
an increased a/c in other dichalcogenides 45] . The a/c 
lattice constant ratio is smaller than in 2H-NbSe2, so 
that pressure drives farther away from the CDW insta- 
bility [45]. This can be fully confirmed, obviously, only by 
measurements under pressure of any property sensitive to 
CDW order, such as the resistance. Nevertheless, if CDW 
re-enters, T C (P) should show some kink or anomaly, and 
a decrease instead of the increase we observe here. Thus, 
our results also show that the re-entrance of CDW is a 
very unlikely possibility. 

It has been postulated that the coexistence of super- 
conductivity with charge density wave is related to the 



neighborhood to a quantum critical point [45J, |48h50| . 

Quantum critical points appear when a second order 
phase transition is driven to zero temperature by modify- 
ing composition or lattice parameters. Quantum fluctu- 
ations diverge at these points and are expected to induce 
emergent exotic properties [5Cj . In the transition metal 
dichalcogenides, quantum critical points may appear hid- 
den below the superconducting or charge ordered states 
[20l ] . Pressure in 2H-NbS2 drives the system farther away 
from such a quantum critical point, which has a marginal 
effect on the critical temperature, and is not associated 
to maxima or dome- like shapes of T c . Our data show 
that maxima in T c can be reached by going farther away 
from CDW instability. Therefore, such dome-like shapes 
in these compounds can be obtained without relation to 
CDW order. 

Finally, let us remark that recent reports discuss syn- 
thesis and characterization of single and few layer sys- 
tems of this and other dichalcogenide compounds. Their 
superconducting properties, as well as the charge den- 
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sity wave, are expected to differ from the bulk 
In these systems, the strain induced in the fabrication 
method is possibly significant and will be of importance 
in the pairing interaction. The pressure dependence of 
the bulk properties should be thus useful to predict and 
understand the modifications found the size reduction 



down to single or few layers 56|. 

In summary, we have presented results on the variation 
of T c and H C 2 as a function of pressure in the dichalco- 
genide 2H-NbS2 . Our data show that the critical temper- 
ature increases smoothly as pressure is increased. This 
behavior is compared to the one already found in 2H- 
NbSe2, where a maximum in T c is found at 10 GPa. 
On the other hand, the upper critical field of 2H-NbS2 
exhibits an intriguing behavior. There is an initial de- 
crease, contrasting with the increase of T c , but above 8.7 
GPa the upper critical field rises again. We provide a 
model to explain this behavior in terms of pressure in- 
duced changes in the Fermi surface. 

The Laboratorio de Bajas Temperaturas is associated 
to the ICMM of the CSIC. This work was supported by 
Spanish MINECO (Consolider-Ingenio CSD2007-00010 
and CSD2009-00013 programs, SAB2009-0057, FIS2011- 
23488 and MAT2011-27470-C02-02), by the Comunidad 
de Madrid through program Nanobiomagnet and by 
NanoSc-COST program. 



[1] P. Monceau, Adv. Phys. 61, 325 (2012). 

[2] T. Yokoya, T. Kiss, A. Chainani, S. Shin, M. Nohara, 

and H. Takagi, Science 294, 2518 (2001). 
[3] E. Boaknin, M. A. Tanatar, J. Paglione, D. Hawthorn, 

F. Ronning, R. Hill, M. Sutherland, L. Taillefer, 

J. Sonier, S. Hayden, et al., Phys. Rev. Lett. 90, 117003 

(2003). 



[4] J. G. Rodrigo and S. Vieira, Physica C 404, 306 (2004). 
[5] T. Kiss, T. Yokoya, A. Chainani, S. Shin, T. Hanaguri, 

M. Nohara, and H. Takagi, Nature Phys. 3, 720 (2007). [31 

[6] I. Guillamon, H. Suderow, F. Guinea, and S. Vieira, 

Phys. Rev. B 77, 134505 (2008). [32 

[7] L. N. Bulaevskri, Sov. Phys. Usp. 19, 836 (1976). 
[8] R. V. Coleman, B. Giambattista, P. K. Hansma, A. John- 
son, W. W. McNairy, and C. G. Slough, Adv. Phys. 37, [33 

559 (1988). 
[9] W. Sacks, D. Roditchev, and J. Klein, Phys. Rev. B 57, [34 

13118 (1998). 
[10] K. Rossnagel, J. Phys.: Condons. Matter 23, 213001 [35 

(2011). [36 

[11] I. Guillamon, H. Suderow, S. Vieira, L. Cario, P. Diener, 

and P. Rodiere, Phys. Rev. Lett. 101, 166407 (2008). [37 

[12] J. Kacmarcik, Z. Pribulova, C. Marcenat, T. Klein, 

P. Rodiere, and L. C. P. Samuely, Phys. Rev. B 82, [38 

014518 (2010). 
[13] A. F. Kusmartseva, B. Sipos, H. Berger, L. Forro, and [39 

E. Tutis, Phys. Rev. Lett. 103, 236401 (2009). 
[14] J. A. Wilson, F. J. DiSalvo, and S. Mahajan, Adv. Phys. [40 

24, 117 (1975). 
[15] J. C. Tsang, M. W. Shafer, and B. L. Crowder, Phys. 

Rev. B 11, 155 (1975). [41 

[16] H. Mutka, Phys. Rev. B 28, 2855 (1983). 
[17] J. L. Vicent, S. J. Hillenius, and H. V. Coleman, Phys. 

Rev. Lett. 44, 892 (1980). [42 

[18] I. Guillamon, H. Suderow, A. Fernandez-Pacheco, 

J. Sese, R. Cordoba, J. M. D. Teresa, M. R. Ibarra, and 

S. Vieira, Nat. Phys. 5, 651 (2009). [43 

[19] E. Coronado, C. Martf-Gastaldo, E. Navarro-Moratalla, 

A. Ribcra, S. J. Blundell, and P. J. Baker, Nat. Chem. 

2, 1031 (2010). [44 

[20] Y. Feng, J. Wang, R. Jaramillo, J. van Wezel, S. Haravi- 

fard, G. Srajer, Y. Liu, Z.-A. Xu, P. B. Littlewood, and [45 

T. F. Rosenbaum, PNAS 109, 7224 (2012). [46 

[21] E. Morosan, H. W. Zandbergen, B. S. Dennis, J. W. G. 

Bos, Y. Onose, T. Klimczuk, A. P. Ramirez, N. P. Ong, [47 

and R. J. Cava, Nat. Phys. 2, 544 (2006). 
[22] B. Sipos, A. F. Kusmartseva, A. Akrap, H. Berger, 

L. Forro, and E. Tutis, Nat. Mater. 7, 960 (2008). [48 

[23] H. Suderow, V. G. Tissen, J. P. Brison, J. L. Martinez, 

S. Vieira, P. Lejay, S. Lee, and S. Tajima, Phys. Rev. B [49 

70, 134518 (2004). [50 

[24] H. Suderow, V. G. Tissen, J. P. Brison, J. L. Martinez, [51 

and S. Vieira, Phys. Rev. Lett. 95, 117006 (2005). [52 

[25] A. Soumyanarayanan, M. M. Yee, Y. He, J. van Wezel, 

D. J. Rahn, K. Rossnagel, E. W. Hudson, M. R. Nor- [53 

man, and J. E. Hoffman, arXiv: 1212.4087 [cond-mat.str- 

cl] (2012). [54 

[26] C. M. Fang, A. R. H. F. Ettema, C. Haas, G. A. Wiegers, 

H. van Leuken, and R. A. de Groot, Phys. Rev. B 52, 

2336 (1995). 
[27] E. Helfand and N. Werthamer, Phys. Rev. 147, 288 [55 

(1966). 
[28] P. C. Hohenberg, N. R. Werthamer, Physical Review [56 

153, 493 (1967). 
[29] L. F. Mattheiss, Physical Review B 1, 373 (1970). 
[30] S. V. Shulga, S. L. Drechsler, G. Fuchs, K. H. Muller, 



K. Winzer, M. Heinecke, and K. Krug, Phys. Rev. Lett. 

80, 1730 (1998). 

T. Dahm and N. Schopohl, Phys. Rev. Lett. 91, 017001 

(2003). 

V. Moshchalkov, M. Menghini, T. Nishio, Q. Chen, 

A. Silhanek, V. Dao, L. Chibotaru, N. Zhigadlo, and 
J. Karpinski, Phys. Rev. Lett. 102, 117001 (2009). 

M. Leroux, P. Rodiere, L. Cario, and T. Klein, Physica 

B Cond Matt 407, 1813 (2012). 

P. Diener, M. Leroux, L. Cario, T. Klein, and P. Rodiere, 

Phys. Rev. B 84, 054531 (2011). 

W. Fisher, M. Sienko, Inorg. Chem. 19, 39 (1980). 

P. D. Horn and Y. M. Gupta, Phys. Rev. B 39, 973 

(1989). 

N. Tateiwa and Y. Haga, J. Phys.: Conf. Series 215, 

012178 (2010). 

B. E. J. Jr., H. H. Shanks, D. K. Finnemore, and B. Mo- 
rosin, Phys. Rev. B 6, 835 (1972). 

P. Molinie, D. Jerome, and A. J. Grant, Phil. Mag. p. 
1091 (1974). 

N. Barisic, I. Smiljanic, P. Popcevic, A. Bilusic, E. Tutis, 
A. Smontara, H. Berger, J. Jacimovic, O. Yuli, and 
L. Forro, Phis. Rev. B 84, 075157 (2011). 
G. Shaw, P. Mandal, S. S. Banerjee, A. Niazi, A. K. 
Rastogi, A. K. Sood, S. Ramakrishnan, and A. K. Grover, 
Phys. Rev. B 85, 174517 (2012). 

M. A. Measson, D. Braithwaite, J. Flouquet, G. Seyfarth, 
J. P. Brison, E. Lhotel, C. Paulsen, H. Sugawara, and 
H. Sato, Phys. Rev. B 70, 064516 (2004). 
L. Glemot, J. Brison, J. Flouquet, A. Buzdin, I. Sheikin, 
D. Jaccard, C. Thessieu, and F. Thomas, Phys. Rev. 
Lett. 82, 169 (1999). 

M. Prohammer and E. Schachinger, Phys. Rev. B 36, 
8353 (1987). 

A. H. C. Neto, Phys. Rev. Lett. 86, 4382 (2001). 
M. D. Johannes, I. I. Mazin, and C. A. Howells, Phys. 
Rev. B 73, 205102 (2006). 

M. Leroux, M. L. Tacon, M. Calandra, L. Cario, M. A. 
Measson, P. Diener, E. Borrissenko, A. Bosak, and 
P. Rodiere, Phys. Rev. B 86, 155125 (2012). 
S. Sachdev, Quantum Phase Transitions (Cambridge 
University Press, Cambridge, 1999). 
P. Coleman and A. J. Schofield, Nature 226, 433 (2005). 
P. Coleman, Phys. Status Solidi B 247, 506 (2010). 
R. Frindt, Phys. Rev. Lett. 28, 299 (1972). 
X. Win, D. Yang, R. Frindt, and J. Irwin, Phys. Rev. B 
44, 3490 (1991). 

A. Castellanos-Gomez, N. Agra'it, and G. Rubio- 
Bollinger, Appl. Phys. Lett. 96, 213116 (2010). 
J. A. Galvis, P. Rodiere, I. Guillamon, M. Osorio, 
G. Rubio-Bollinger, J. Rodrigo, L. Cario, E. Navarro- 
Moratalla, E. Coronado, S. Vieira, et al., ArXiV 
1210.2659 (2012). 

Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, 
and M. S. Strano, Nature Nanotechnology 7, 699 (2012). 
M. Calandra, I. I. Mazin, F. Mauri, Physical Review B 
80, 241108(R) (2009). 



